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Understanding the electronic properties of dopants near an interface is a critical challenge for 
nano-scale devices. We have determined the effect of dielectric mismatch and quantum confinement 
on the ionization energy of individual acceptors beneath a hydrogen passivated silicon (100) surface. 
Whilst dielectric mismatch between the vacuum and the silicon at the interface results in an image 
charge which enhances the binding energy of sub-surface acceptors, quantum confinement is shown 
to reduce the binding energy. Using scanning tunneling spectroscopy we measure resonant transport 
through the localized states of individual acceptors. Thermal broadening of the conductance peaks 
provides a direct measure for the absolute energy scale. Our data unambiguously demonstrates that 
these two independent effects compete with the result that the ionization energy is less than 5 meV 
lower than the bulk value for acceptors less than a Bohr radius from the interface. 



The operation of semiconductor devices is based on the 
possibility to locally change the electron properties of the 
host material by means of doping. As device dimensions 
continue to decrease, the surface- to- volume ratio of active 
channels increases and the effect of the semiconductor- 
insulator interface on local doping starts to dominate de- 
vice properties [1 j. Previous studies have suggested that 
dielectric mismatch at the semiconductor-insulator inter- 
face leads to an increase in ionization energy of dopants 
near the interface [2j [3] . In silicon nanowires this leads 
to doping deactivation and consequently an increase of 
resistivity with decreasing diameter [2j [4]. However, re- 
cent transport spectroscopy experiments on single arsenic 
donors in gated nanowires did not report an appreciable 
increase in ionization energy [5 j. These results appear 
contradictory. 

The ionization energies of shallow donor and acceptor 
impurities are qualitatively described by effective mass 
theory [6 , which works especially well for light impurity 
atoms such as Li and B. Since the Coulomb potential 
is strongly screened due to the polarization of the semi- 
conductor, the ionization energy of dopant impurities is 
only in the order of tens of meV. This simple picture 
breaks down in the presence of an interface. Dielectric 
mismatch between the semiconductor material and its 
surroundings is predicted to enhance the ionization en- 
ergy [4 j. Moreover, for nanowires it is well known that 
when the thickness of the nanowire approaches the Bohr 
radius of the impurity the ionization energy increases due 
to quantum confinement [7J[8]. However, in the case of a 
half-space, i.e. a flat interface, effective mass theory pre- 
dicts a decrease in the ionization energy due to quantum 
confinement [9j [10] . As a result of these two competing 
effects, dielectric mismatch and quantum confinement, 
the ionization energy of dopant near a flat interface is 
expected to be bulk-like [9] [10] . 

Here, we use low-temperature (4.2 K) scanning 



tunneling spectroscopy (STS) to directly measure the 
ionization energy of boron acceptors beneath the hydro- 
gen terminated Si(100) surface (Na ~ 8 x 10 18 cm -3 ). 
Previous studies on GaAs(llO) [HHIS], InAs(llO) [E] 
and ZnO(0001) [17] surfaces have proven that STS is 
a powerful tool to study sub-surface impurities. Using 
scanning tunneling microscopy the surroundings of each 
individual dopant atom can be imaged and therefore 
any effect of dopant clustering or interface disorder, 
such as charge traps, on the ionization energy of the 
acceptors may be excluded. By analyzing the line 
shape of differential conductance within well known 
single electron spectroscopy formalisms, we conclude 
that transport is thermally broadened. Consequently, 
we can use the electron temperature of the sample 
as a reference to calibrate the coupling between the 
applied bias voltage and the potential landscape at 
the semiconductor- vacuum interface. In such a way we 
are able to obtain a direct measure for the acceptor 
ionization energy. Moreover, STS allows us to determine 
the distance of individual acceptors to the interface by 
measuring the spectral shift of the valence band due to 
the negatively charged acceptor nucleus. Importantly, 
from the thermally broadened single-electron transport 
through the localized acceptor state, in conjunction with 
the spectral shift of the valence band due to the ionized 
nucleus, the dopant depth can be directly correlated to 
its ionization energy. 

The schematic energy diagrams in Fig. [I] illustrate 
three different transport regimes in which we study the 
sub-surface boron acceptors: (i) charge sensing in the 
valence band, (ii) resonant tunneling through the local- 
ized acceptor state and (iii) charge sensing in the con- 
duction band. When a negative sample bias voltage is 
applied, the presence of a sub-surface acceptor results in 
an increase in the direct tunneling into the valence band 
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FIG. 1. Schematic energy diagram of the tip- vacuum-sample 
tunnel junction, (i) When eV < Ey < cVfb electrons tun- 
nel from the valence band to the tip and the presence of a 
sub-surface acceptor is observed as a protrusion in the STM 
topography, (ii) In case cVfb < eV ~ Ec direct tunneling 
via the localized acceptor state leads to a protrusion in the to- 
pography, (iii) Supression of the local density of states by the 
acceptor potential reduces the tunnel current when Ec < eV, 
resulting in a depression in the topography. 



due to an increase in the local density of states (LDOS) 
caused by the negatively charged nucleus of the acceptor. 
The increase in tunnel current leads to a height increase 
in the STM topography, as is shown in figure Qi). Like- 
wise, when a positive voltage is applied, suppression of 
the local density of states by the acceptor potential leads 
to a decrease in the direct tunnel current into the con- 
duction band and resulting in a dip in the STM topog- 
raphy [16] (Fig. [TJiii)). However, when a positive voltage 
is applied such that the Fermi energy of the tip is close 
to the conduction band edge, the transport is no longer 
dominated by direct tunneling into the conduction band 
but by resonant transport through the localized acceptor 
state (Fig[jjii)). 

Since the Fermi energy of the heavily B-doped sample 
is pinned at the bulk acceptor level E^ lk , the voltage at 
which the localized acceptor state is equal to the Fermi 
level, Vonseti with respect to the fiat-band voltage, Vfb, 
is a direct measure for the difference in ionization energy, 
AE = Ea — E^ lk : of the sub-surface acceptor with re- 
spect to the bulk ionization energy. The bias dependence 
of the energy level of the localized acceptor state Ea is 
described by the lever arm a = e (IEa/cIV [HI [19]. In 



this Letter we present a direct measurement of the lever 
arm and measure the shift AE = —ae(V onset — Vfb) 
by studying transport through individual sub-surface 
acceptors where we directly determine: (i) the potential 
due to the negatively charged nucleus, (ii) the flat-band 
voltage and depth of the acceptors from direct tunneling 
to/from the conductance/ valence band and (iii) the on- 
set voltage and lever arm from single-electron transport 
through the localized acceptor state. 

Figure [g|a) shows the normalised conductance Gn = 
(dI/dV)/(I/V) measured away from any sub-surface ac- 
ceptor. The flat-band voltage Vfb was extracted by 
comparing the 4.2 K bandgap of Si(001):H with voltages 
Vy and Vc for tunneling into the valence and conduc- 
tance band edges (dashed lines in Fig. [2|a) and (b)). A 
first approximation of Vfb is made by taking the tip 
induced band bending (TIBB) to be linear (dotted line 
in Fig. |2jb)). This constraint is subsequently relaxed in 
order to account for screening (solid line in Fig. [2jb)). 
Finally, the flat-band condition was independently mea- 
sured from the apparent barrier height. 

The onset voltage for tunneling from the valence band 
Vy and tunneling into the conduction band Vc is deter- 
mined by finding the voltage axis intercept (i.e., Gn = 0) 
of the linear extrapolation of the Gn — V curve at its max- 
imum slope point (dotted lines in Fig. [2|a)) [20]. The 
potential (j)siV) at the interface as a function of the ap- 
plied sample voltages is obtained by from the flat-band 
energies E F - Ey = E b / lk - Ey = 0.045 eV [21 and 
E g = E c - E v = 1.17 eV at T = 4.2 K [22] away from 
the acceptor. Assuming a linear relationship between the 
applied voltage and the potential at the interface 



dfc = e(V c - Vy) - E g 
dV Vc-Vv 



(1) 



yields an approximated flat-band voltage Vfb = 0.5 ± 
0.1 V which is inferred from the condition for tunneling 
from the valence band 4>siVy) = eVy 



d<j>s 
dV 



(Vy - Vfb) = eVy. 



(2) 



When the flat-band voltage lies within the band gap the 
d(j>s/dV can not be assumed linear as it is well known 
that for V > Vfb accumulated carriers at the surface will 
screen the electric field from the tip more strongly. In or- 
der to correct for this effect we use the <j>s{Y) calculated 
using the method of Feenstra [23 (Fig. |5Jb)) and deter- 
mine a corrected flat-band voltage Vfb — 0.38 ± 0.1 V. 
This flat-band voltage is smaller than expected from the 
difference between the bulk workfunction $ = 4.55 eV 
[24] of tungsten and sample electron affinity x — 4.05 eV 
[25] . The measured flat-band voltage corresponds to a 
tip workfunction $> tip = 4.9 ± 0.1 eV (larger values for 
the tip workfunction have been previously reported) . Fi- 
nally, we independently confirmed the value of the tip 
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FIG. 2. (a) Normalized conductance Gn measured on a hy- 
drogen terminated Si(100):H surface away from (solid line) 
and above (dashed line) an acceptor. The top of the valence 
band Vy and bottom of the conduction band Vc are deter- 
mined by the slopes oi Gn away from the acceptor (dotted 
lines), (b) Surface potential <j)s assuming linear TIBB (dot- 
ted line) and non-linear TIBB (solid line) calculated using the 
code developed by Feenstra [23 . 



workfunction by measuring the apparent barrier height 
P~4j|26] ( see section 1 of the Supplementary Material). 

The distance of the sub-surface acceptors to the inter- 
face is measured from the spectral shift of the valence 
band edge. Figure [3^ a) shows the differential conduc- 
tance (dl/dV) map, as a function of position and sample 
bias voltage V, measured simultaneously with the empty 
state topography at V = 2 V shown in Fig. |3ja). At the 
acceptor site the dl/dV map clearly shows an upward 
spectral shift of the valence band states due to the buried 
acceptor. The spatially resolved shift of the valence-band 
edge AEy as shown in Fig. [3jb) is defined from the slope 
of Gat as before. The first-order perturbation to the bind- 
ing energy of the valence band states at the interface, ip 3 , 
and thus the shift of the valence-band edge, AEy, due to 
the potential of the negatively charged acceptor nucleus 
at position ro can be estimated as 

AE v = {il> a \U A {r,r )\il>s), (3) 

where Ua(t, ro) is the acceptor potential at position r. 
In the case of a classical half-space, i.e. in the absence 
of screening by free carriers in the tip or the substrate, 
the approximation (^s|?7a(^j r o)\^s) ~ ^(^ro), where 
s is the lateral separation with respect to the acceptor 
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FIG. 3. (a) dl/dV map recorded simultaneously with empty- 
state topography (+2.2 V, 300 pA). Tunneling from the va- 
lence band to the tip is indicated by VB, tunneling from the 
tip to the conduction band is indicated by CB. The dl/dV 
map is cut at the acceptor site to show the upward shift of 
the valence band states, (b) Schematic potential landscape 
due to a negatively charged nucleus below the sample surface 
and its image charge in the vacuum. The solid line indicates 
the potential along the vacuum-semiconductor interface, (c) 
Shift in the valence band maximum Ey as function of lat- 
eral tip position, (d) Azimuthal averaged shift in the valence 
band maximum Ey (filled squares) as function lateral tip sep- 
aration s from the acceptor as indicated in (c). The dopant 
depth is determined by fitting Ey to a bare Coulomb poten- 
tial (solid line), taking into account the image charge due to 
dielectric mismatch as illustrated in (b). 

nucleus along the interface as shown in Fig. j3^b), yields 

e 2 1 

AE V « . (4) 

47re e e// ^s 2 + d 2 V 7 

where the modified dielectric constant e e ff and depth 
d can be independently determine from a fit [I5j [27] . 
The modified dielectric constant at the interface 
e e ff = (e v + esi)/2 is due to the mismatch between the 
dielectric constants e v and es% = 11.4 [28] of the vacuum 
and silicon, respectively, which leads to a single image 
charge at — d as shown in Fig. |3jc). Figure |3jd) shows 
measured (filled squares) and fitted (solid line) spectral 
shift of the valence-band edge AEy as function of lateral 
tip separation s from the dopant. Following references 
P~5j [27] we fit the spectral shift of the valence-band 
edge as a function of position to equation [4] (Fig. (3^d)) 
and independently determine the depth d of individual 
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acceptors and the modified dielectric constant e e ff- 
Importantly, any screening by either by carriers in 
the tip or in the substrate would reduce the effect of 
dielectric mismatch and thus lead to a deviation of the 
modified dielectric constant from e e // = (e v + esi)/2. 
The obtained modified dielectric constant for all five 
measured acceptors agree within experimental error 
with the expected value e e // = 6 following the classical 
half-space approach and experimental values that have 
previously been reported for STM experiments p~5j [27] . 

We use the thermal broadening of the conductance 
peaks to calibrate the lever arm a for individual accep- 
tors. The shift of the acceptor energy levels AE due 
to the applied bias voltage V depends on the screening 
of the electric field in the semiconductor and the over- 
lap between the acceptor wavefunction and the TIBB 
(Fig. ga)). When screening is strong, the majority of 
the applied voltage will be dropped over the vacuum sep- 
arating the tip and the sample and the potential in the 
semiconductor will only be slightly modified by the ap- 
plied sample voltage. Furthermore, the lever arm will 
be smaller for acceptors deeper below the surface, as the 
electric field is more strongly screened further from the 
interface. 

In previous studies the tip induced band bending has 
been estimated by solving the one- or three-dimensional 
Poisson equation p~5J [19] [23j [26]. However, the calcu- 
lated potential in the semiconductor will strongly depend 
on experimental parameters, e.g. tip-sample separation, 
tip shape and carrier concentration. Although in some 
cases these experimental parameters can by fitted from 
the data [19], they are typically not known. Moreover, 
the estimate of the lever arm from calculating the over- 
lap between the TIBB and the impurity wavefunction is 
very sensitive to exact size of the wavefunction, which is 
well known to differ from the bulk wavefunction [I8j [19] . 
Rather than trying to estimate the lever arm by solving 
the Poisson equation, we fit the conductance peak to a 
thermally broadened Lorentzian [29], as shown Fig.gb). 
This line shape describes resonant tunneling via a single 
localized state into a continuum of states that are occu- 
pied according to the Fermi-Dirac distribution [30j [31] . 
Other broadening mechanisms, e.g. lifetime broadening 
and charge noise, can be excluded as dominant sources of 
conductance peak broadening since they do not have the 
appropriate line shape of the conductance peaks and the 
fact that the width of the peaks does not increase with 
sampling time (see section 2 of the Supplementary Mate- 
rials). By allowing a to vary linear with V in our fit we 
observe that the measured values for a do not depend on 
the bias voltage V. Consequently, we can conclude Stark 
shifts of the localized states due to the electric field is 
negligible in our measurement geometry. For consistency 
we compare our method with the previously described 
method [I5J [T9l E2 [26] . We have calculated V onse t and 
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FIG. 4. (a) Schematic energy diagram. Local tip induced 
band bending brings the acceptor level Ea into resonance 
with the Fermi energy. The voltage Vonset at which this oc- 
curs depends on the acceptor depth, ionization energy and the 
screening length, (b) The lever arm a and the onset voltage 
Vonset are determined fitting the two conductance peaks in the 
bandgap to the sum of two thermally broadened Lorentzian 
line shapes, (c) Measured (squares) and calculated (line) on- 
set voltage for resonant tunnelling as a function of acceptor 
depth, (d) Measured (squares) and calculated (line) lever 
arm as a function of acceptor depth, (e) The shift in ioni- 
sation energy AE with respect to the bulk ionisation energy 
inferred from the measured onset voltage and the lever arm 
(squares). The onset voltage and lever arm are calculated 
for AE — —2.5 meV, line in (c), and an overlap between 
the acceptor wavefunction and the tip induced band bending 
AE = -0.06 x TIBB, where TIBB is the tip induced band 
bending at the acceptor site as illustrated in (a). 



a as a function of depth for AE = 2.5 meV where we 
use the overlap between the impurity wavefunction and 
the TIBB as a free parameter to match the empirical re- 
sult. We find that for AE = -0.06 x TIBB the Poisson 
model matches our results. Importantly, the 0.06 over- 
lap is a factor ^4 smaller the 0.27 overlap predicted in 
previous studies [T3l [19] . The predicted uncertainty for 
the Poisson model is of the order of a factor 4 [15] . 

Figure gc) shows the onset voltage V onset as a 
function of dopant depth for five different acceptors. 
The measured lever arm a is shown if Fig. gd). The 
shift in ionization energy AE = —ae(V onse t — Vfb) 
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is smaller than zero for all measured acceptors, as 
shown in Fig. [4^e). Since the measured value of e e ff 
corresponds to the expected value following the classical 
half-space approach we can conclude that two opposing 
effects influence the ionization energy of near-interface 
dopant atoms: (i) dielectric mismatch; the dopant 
potential, which is screened by charge polarization 
in the semiconductor, becomes more attractive when 
its environment becomes less polarizable, i.e. has a 
lower dielectric constant, which leads to an increase 
of the ionization energy; (ii) quantum confinement; 
exclusion of the dopant wavefunction from the region 
outside the semiconductor results in a decrease of 
the ionization energy [9j [10]. The observed bulk-like 
ionization energies for acceptors less than an effective 
Bohr radius from the interface are strong evidence 
that the effect of dielectric mismatch at the interface 
is mitigated by quantum confinement. Our transport 
data unambiguously demonstrates that acceptors within 
an effective Bohr radius from the interface of a silicon 
half-space geometry are not deactivated. We would like 
to point out that atomistic differences from the bulk 
may lead to an alteration of the binding energy such as 
the enhancement observed in the same geometry in Ref. 
[18], but this effect is unrelated to dielectric mismatch. 

In conclusion, we have determined the ionization 
energy of individual sub-surface acceptors below the 
Si(100):H surface by means of low-temperature scanning 
tunneling spectroscopy. Calibration of the local lever arm 
using the thermal broadening of the conductance peaks 
removes the necessity of modeling the tip induced electro- 
static potential below the semiconductor surface, provid- 
ing a direct measure for the ionization energy. We believe 
that this method provides a valuable parallel between ex- 
tensively studied electron transport in mesoscopic devices 
and scanning tunneling spectroscopy. Moreover, this ex- 
periment demonstrates bulk-like ionization energies for 
acceptors less than a Bohr radius away from the inter- 
face. The fact that sub-surface acceptors are not deacti- 
vated is of great importance for the doping of nanoscale 
devices. 
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